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A  Semidefinite  Programming  Approach  to  Source 
Localization  in  Wireless  Sensor  Networks 

Chen  Meng,  Student  Member,  IEEE ,  Zhi  Ding,  Fellow,  IEEE ,  and  Soura  Dasgupta,  Fellow,  IEEE 


Abstract — We  propose  a  novel  approach  to  the  source  local¬ 
ization  and  tracking  problem  in  wireless  sensor  networks.  By 
applying  minimax  approximation  and  semidefinite  relaxation,  we 
transform  the  traditionally  nonlinear  and  nonconvex  problem 
into  convex  optimization  problems  for  two  different  source  local¬ 
ization  models  involving  measured  distance  and  received  signal 
strength.  Based  on  the  problem  transformation,  we  develop  a  fast 
low-complexity  semidefinite  programming  (SDP)  algorithm  for 
two  different  source  localization  models.  Our  algorithm  can  either 
be  used  to  estimate  the  source  location  or  be  used  to  initialize  the 
original  nonconvex  maximum  likelihood  algorithm. 

Index  Terms — Maximum  likelihood  estimation,  semidefinite 
programming,  source  localization,  wireless  sensor  network. 


I.  Introduction 

THE  problem  of  locating  a  signal  source  using  a  wireless 
network  of  sensors  has  been  addressed  in  several  past 
papers  [1],  [2]  and  references  therein.  One  solution  to  this 
problem  applies  weighted  least  squares  (WLS)  estimation 
[3]— [5]  that  can  be  quite  effective.  However,  the  WLS  approach 
operates  with  a  cost  function  that  is  nonconvex  and  subject 
to  local  minima,  thus  requiring  either  smart  initialization  or 
special  conditions  to  ensure  global  convergence.  To  remedy 
the  local  convergence  problem,  an  algorithm  based  on  the 
concept  of  projection  onto  convex  sets  (POCS)  has  recently 
been  proposed  in  [6]  and  [7].  The  POCS  algorithm  is  very 
effective  in  overcoming  the  local  convergence  problem  when 
the  source  resides  within  the  convex  hull  of  the  sensor  nodes. 
Nevertheless,  the  POCS  algorithm  does  suffer  from  poor  per¬ 
formance  when  trying  to  locate  a  source  node  not  inside  the 
convex  hull.  In  this  work,  we  present  a  new  approach  to  tackle 
the  local  convergence  problem  that  the  WLS  approach  faces 
and  the  problem  with  source  outside  the  convex  hull  that  the 
POCS  algorithm  faces.  By  transforming  the  source  localization 
problem  into  a  convex  optimization  problem  using  minimax 
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approximation  and  semidefinite  relaxation,  we  can  find  the 
global  minimum  of  the  relaxed  problem  and  do  not  need  to 
impose  coverage  conditions. 

The  application  of  semidefinite  programming  (SDP)  in  sensor 
localization  with  noisy  distance  measurement  and  angle  of  ar¬ 
rival  has  been  studied  by  previous  researchers,  e.g.,  [8]  and  ref¬ 
erences  therein.  Nevertheless,  instead  of  minimizing  the  l1  mea¬ 
surement  error,  this  work  derives  the  SDP  from  the  maximum 
likelihood  estimation  and  the  minimax  approximation.  Another 
important  feature  of  our  algorithm  is  that  it  can  be  applied  to  two 
sensor  measurement  models.  In  other  words,  the  SDP  can  be  ap¬ 
plied  to  the  problem  of  measured  distance  with  additive  white 
Gaussian  noise  and  the  problem  of  received  signal  strength  with 
lognormal  fading. 

Linally,  the  newly  proposed  SDP  algorithm  can  be  efficiently 
implemented  via  a  standard  numerical  optimization  toolbox 
such  as  SeDuMi  [9].  The  implementation  also  has  very  fast 
convergence.  Indeed,  the  computation  complexity  is  polyno¬ 
mial-time  [10],  with  the  worst-case  scenario  of  0(m6  5),  where 
m  is  the  dimension  of  space  in  which  sensors  are  placed.  This 
feature  makes  it  appealing  for  resource-limited  applications 
and  for  tracking  targets  in  real  time. 

II.  Problem  Statement 

We  now  describe  our  problem  formulation  and  basic  assump¬ 
tions.  The  known  sensor  positions  are  a  set  of  m— dimensional 
vectors  xi, . . . , Xjy-  The  signal  source  location  is  an  unknown 
vector  y.  An  source  location  estimate  y  is  determined  from 
sensor  measurements.  In  this  letter,  we  only  consider  single 
measurement  from  each  sensor  node.  Nevertheless,  our  results 
are  readily  applicable  to  cases  involving  multiple  measure¬ 
ments.  We  consider  the  following  two  measurement  models. 
Note  that  the  major  difference  between  the  two  models  is  the 
different  assumption  of  noise,  and  the  scenarios  can  be  found 
in  the  referred  literature. 

A.  Distance  Measurement  Model 

In  this  common  model  [5],  [11],  the  distance  di  from  source 
to  each  sensor  node  is  measured,  i.e., 

df  =  ||x;  -y||2  +ra*,  i  =  l,...,N  (1) 

where  the  noise  nt  is  i.i.d.  Gaussian  with  zero  mean  and  variance 
a2.  The  conditional  probability  density  is 

i  E i'2)2 

p(4,4,...,4r  |y)  =  (27r(j2)Jv/2e 
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and  the  maximum  likelihood  estimate  (M  L E)  of  y  is 

N 

y  =  argmin^]  (d?  -  ||x;  -y||2)2.  (2) 

y  i- 1 

This  model  is  less  realistic  in  that  nt  cannot  be  truly  Gaussian 
as  dj  is  nonnegative. 

B.  Signal  Strength  Measurement  Model 

If  the  source  at  y  emits  a  signal  of  power  P,  the  sensor  at  x; 
can  receive  a  signal  power  s; 


By  denoting  Y  =  yyT  and  applying  semidefinite  relaxation 
[10],  [14],  this  problem  can  be  reformulated  into 

min  t 

s.t.  —  t  <  Trace (x.,Y)  —  d2  <t 
i  =  1,2,...,  N 

Y  y  0,  Y(m  +  l,m  +  1)  =  1  (7) 

where  Y  y  0  denotes  (symmetric)  positive  semidefinite.  Equa¬ 
tion  (7)  is  a  convex  optimization  problem.  Its  global  optimal  so¬ 
lution  can  be  found  using  modern  5DP  solvers  such  as  SeDuM  i 
[9]  that  applies  the  interior  point  method. 


Si  =  P||xi  -  y ||  p  (3) 

where  (3  is  the  path  loss  coefficient  [12].  Under  lognormal 
fading,  the  received  signal  strength  (RSS)  in  dB  at  a  sensor  can 
be  modeled  as  [7],  [12],  [13] 


10  log  Si  =  10  log  P  —  10/31og(||xj  -  y  || )  +  n*  (4) 


B.  Semidefinite  Relaxation  for  Received  Signal  Strength 

We  first  reformulate  the  localization  problem  of  (5)  into  a 
minimax  problem 


y  =  argmin  max 

J  &  J=r|. v 


=  argmin  max 
i=l,2...,N 


(8) 


where  the  noise  m  is  assumed  i.i.d.  Gaussian  with  zero  mean 
and  variance  a2.  The  conditional  probability  density  is 

p(lnsi,lns2,...,lnsjv|y) 

=  _ p  2(ln  10/10-cr)2 

(27r(ln  10/10  •  a)2)N/2 
and  the  M  LE  of  y  can  be  obtained  as 

y  =  arg mm  E  (ln  ~  ^  ( |[x.  Ty y/i ) )  •  (5) 


We  can  show  that  this  new  optimization  problem  is  less  sen¬ 
sitive  to  the  path  loss  parameter  6.  To  do  so,  first  note  that  from 
(4),  we  have 


— 


p  .  ion^/10 
l|x.t-y||^ 


The  minimum  of  the  cost  function  in  (8)  assumes  the  value 


A  y)  = 


In  A 


(*i\W  I, 

(p)  l|x--; 


In  10  -  2  In  ||x*  -  y||  +  2 In  ||x*  -  y|| 


III.  Our  New  Efficient  Algorithms 

A.  Semidefinite  Relaxation  With  Distance  Measurement 

The  optimization  problem  of  (5)  involves  high-order  power 
of  y  and  is  highly  nonlinear  and  nonconvex.  To  simplify,  we 
propose  a  minimax  approximation 

y*  =  argmin  ._max  | d2  -  ||xj  -  y||2| .  (6) 

This  approximation  is  supported  by  the  so-called  equivalence 
between  l2  and  Z°°  vector  norms. 

To  proceed,  define  an  (m  +  1)  x  1  vector  y  =  [yT  l]T  and 

_  _  I  —  Xj 

T'  _ v  /  „T  v . 

Then  (6)  can  be  written  in  an  equivalent  form 

min  t 

s.t.  y (to  +  1)  =  1,  -t<  yTXif  ~  A  <  t,  i  =  1, ...  ,1V 


for  a  certain  i.  The  estimate  y  demonstrates  insensitivity  to  the 
path  loss  parameter  in  that  for  low  noise  m,  small  perturbation 
of  fit  has  little  effect  on  the  estimate  of  y.  Particularly  when  in 
noise-free  cases,  the  minimum  cost  is  independent  of  /?.  This 
is  not  the  case  for  algorithms  such  as  POC5  in  which  8  affects 
the  radius  of  a  convex  set.  Empirically,  (3  only  represents  the 
average  path-loss  effect.  In  a  practical  environment,  different 
sensors  may  have  significantly  different  values  of  (3,  which  is 
a  well-recognized  drawback  of  applying  the  RSS  signal  model 
when  the  path  loss  parameters  vary.  Hence,  algorithms  that  are 
insensitive  to  j3  are  generally  desirable. 

To  find  an  efficient  solution,  first  denote 


We  transform  (8)  into  its  equivalent  form 

min  t 

s.t.  y(m  +  1)  =  1 

-t  <  ln(gj||xi  -  y ||2)  <  t,  i  =  1, 2, . . . ,  N. 
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Using  the  same  relaxation  technique  as  in  the  previous  section, 
we  have  the  new  problem  as 

min  t 

s.t.  —  t  <  ln(</;  •  Trace(xiY))  <  t,  i  =  1, 2, . . . ,  N  (9) 
Y  y  0,  Y(to  +  1,  m  +  1)  =  1. 

To  make  first  constraint  in  (9)  more  tractable,  additional  approx¬ 
imations  can  be  applied.  First  observe  that  In  a;  <  x  -  l.Thus, 
a  necessary  condition  for 

-t  <  ln(®Trace(xiY)) 

is 


-t  <  ®Trace(xiY)  -  1. 

We  therefore  have  two  necessary  conditions 

1  -  t  <  t/,Trace{x,  Y} 

0  <  TracejxiY}. 

Because  both  Y  and  Xi  are  symmetric  positive  semidefinite 


Tracefx^Y}  >  0 

is  guaranteed  and  becomes  redundant.  Now  consider  the  right- 
hand  side  of  the  first  constraint 


ln(® Trace{xi Y })  <  t. 

We  have  a  sufficient  condition 


ln^Tracejx^Y})  <  ^Trace-f^Y}  —  1  <  t. 

To  summarize,  we  transform  (9)  into  an  approximation 

min  t 

s.t.  —  t  <  qi  ■  Trace(xiY)  —  1  <t,  i  =  1, 2, . . . ,  N 

y  y  o,  Y(m  + 1,  to  + 1)  =  l. 

(10) 

The  optimization  problem  in  (10)  can  now  be  solved  using  stan¬ 
dard  SDP  tools.  Note  that,  in  general,  semidefinite  relaxation  in¬ 
creases  the  problem  size,  and  some  postprocessing  techniques 
are  required  to  convert  the  SD  relaxation  solution  Y  into  an  ap- 
proxi  mate  sol  uti  on  of  the  ori  gi  nal  opti  mizati  on  probl  em  i  n  terms 
of  y.  Such  standard  techniques  have  been  routinely  applied  in 
signal  processing  and  communications  [14],  [15]. 

IV.  Simulations 

We  provide  two  examples  for  performance  comparison. 
Example  1:  Consider  the  distance  measurement  model.  Sen¬ 
sors  are  located  in  2-D  at  Xi  =  [2,2]T,x2  =  [2,6]T,x3  = 

[6. 2] T.  The  source  is  first  placed  at  y  =  [0,0]T  (outside  the 
convex  hull  formed  by  the  sensors)  and  later  placed  at  y  = 

[3. 3] t  (inside  the  convex  hull).  Note  that  when  the  source  is 
outside  the  convex  hull,  the  POCS  performance  is  generally  bad 


Fig.  2.  The  3-D  results  under  the  received  signal  strength  model. 


because  the  intersection  of  convex  sets  is  quite  large,  giving 
rise  to  many  potential  solutions.  Also  note  that  in  this  case,  a 
spurious  stationary  pointy'  =  [6,6]T  is  locally  stable,  when 
solving  the  maximum  likelihood  estimation  using  gradient  de¬ 
scent  methods  [5].  In  Fig.  1,  we  compare  the  estimation  results 
of  POCS  and  SDP  for  the  two  different  source  locations.  The  re¬ 
sults  are  given  as  the  mean-squared  error  of  the  estimate  versus 
standard  deviation  of  the  noise,  averaged  over  3000  M  onte  Carlo 
tests.  Good  estimation  is  demonstrated  consistently  for  the  SDP 
method,  while  POCS  fails  to  locate  the  source  consistently. 

Example  2:  We  now  consider  the  RSS  model.  We  set 
p  =  1000,/?  =  3,  and  distribute  five  sensors  in  3-D  at 
xi  =  [1,0,0]t,x2  =  [0, 2,  0]t,x3  =  [— 2,  — 1,0]t,x4  = 
[0,0,2]T,x5  =  [0,0, -l]T.  We  first  put  the  true  source  at 
y  =  [-1,1,  l]T  (outside  the  convex  hull  of  sensors)  and  then  at 
y  =  [o,0, 0]T  (inside  the  convex  hull).  In  both  cases,  the  SDP 
algorithm  provides  consistently  good  estimates  while  POCS 
works  wel  I  for  the  second  source  but  fai  I s  to  I  ocate  the  first  one. 
The  results  are  illustrated  in  Fig.  2. 
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V.  Conclusion 

We  propose  a  novel  approach  to  the  source  localization 
problem  in  wireless  sensor  networks.  Our  algorithm  relies  on 
the  concept  of  minimax  approximation  to  the  optimal  max¬ 
imum  likelihood  estimation,  and  it  takes  the  form  of  an  efficient 
semidefinite  programming.  We  study  two  different  measure¬ 
ment  models  for  source  localization  and  present  well-posed 
solutions  in  both  cases.  We  present  simulation  examples  that 
demonstrate  the  good  performance  and  consistency  of  the  new 
SDP  algorithms.  Our  algorithm  is  also  less  sensitive  to  inac¬ 
curacy  in  the  path  loss  factor  of  the  RSS  model.  The  proposed 
SDP  estimate  can  serve  either  as  a  standalone  solution  or  as 
an  initialization  to  the  more  complex  maximum  likelihood 
algorithms. 
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